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Methods for computing the maximum and minimum of a polynomial with real coefficients in the 
interval [0, 1] are described, and certain bounds are given. 
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Introduction 

Let p(x) = ao+ d\X -f . . . + a n x n have real coefficients. If / is the interval [0, 1] and 

A/=max p(jt), ro = min p(x), 

xel .re I 



we wish to study some methods of approximating m and M relatively easily. Apart from their in- 
trinsic interest, such methods would seem to have application in computations using interval 
arithmetic in which a basic operation is the determination of the range of a rational function with a 
given interval as domain. The restriction of our discussion to the interval / entails no loss of gen- 
erality, since any interval [a, b] can be mapped onto / by a linear function of x, and such changes 
of variable leave the set of polynomials of degree at most n invariant. 

Some estimates of m and M were given by Cargo and Shisha. 1 They first observe that if n ^ 1 
and ./ = {!, . . .,*}, 



mm 



HHi*. 



M- 



I 



m £ p Wi^ M 



(1) 



Then they note that p has a representation in Bernstein form, namely, 






(2) 



which leads immediately to the bounds 



mm bj^ m : 



. M ^ max bj, 

jeN 



(3) 



where /V={0, . . ., n}. 



*An invited paper. 
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1 Cargo, G. T., and Shisha, O., The Bernstein form of a polynomial, J. Res. Nat. Bur. Stand. (U.S.), 70B (Math. Sci.) No. 1, 79-81 (Jan-Mar. 1966). 
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The numbers bo, . . ., b n are determined by 



6 J = 2 ai 7^ J=°> •••.». (4) 



Our work is intended to be a sequel to the paper of Cargo and Shisha (see footnote 1). In 
section 1, we give a sharper result than (1). Section 2 is devoted to a generalization of the repre- 
sentation (2) and a subsequent improvement of (3), including estimates of the precision of the 
bounds obtained for m and M 9 which are due, essentially, to S. Bernstein. 2 The paper concludes 
with section 3, containing some results in the case that the coefficients a , . . . , a n are complex 
numbers. 

1. Estimates Using Function Values 

Theorem 1: Suppose = t <ti< . . . <t k = l, and d k = max (t j+ i — tj), j = 0, . . ., k — 1, 
then, if K={0, . . ., k} 

d 2 d 2 

min p(ti)--^max |p"(x)| ^ m ^ M ^ max p(ti) + -^max |p"(x)|. (5) 

ieK O x ei ieK O x ei 

Proof: Suppose p(ij)=M and 

to-fM*<-fl.i-o, . . ., k, 



d k 



so that 

According to Taylor's formula, we have 

P(t 3 ) =p(€) + (h-€)p'({) + itj ~ 2 i)2 p"(v) (6) 

where rje/. If f = 0, 1, the right-most inequality in (5) is trivially true. If < £ < 1 then p'(f ) =0 
and (6) implies that 



c) 2 d 2 k 

M ^ max p(tj) H max \p" (77) | ^ max p(tj) +— max \p"(j))\. 



(tj-tV _ 

jeK 2 Vel jeK 8 lei 



An entirely analogous argument establishes the lower bound on m in (5). 
COROLLARY: 7/tj = i/k then d k =l/k, and since 

max|p>)|^2(j-l)j| aj | 
(5) now becomes, 

^K np (k)"8fe 2 G-l)llail*m*M* max P^+^S (j-l»|*,|. (7) 

j=0 j=l 



2 Bernstein, S. N., Collected Works, Vol. I, Translation: AEC-tr-3460, pp. 68-73. 
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The improvement of the error estimate in (7) over that in (1) is evident. 

If we consider p(x) = x — x 2 then M = 1/4, which is precisely the upper bound given by (7) with 
k=\. 

2. The Bernstein Form 

Suppose 5 = 0, 1, . . . , k then we have 
Lemma 1: 



k 



Proof: 






x'=x'(x+ (1 -*))*-= 2) (* • S )x s+i {l-x 



)*- 



■|(J:;)«M-.)«-|gj(*)-a 



-*)*- 



Thus, the functions 



(*)*»(! -*)*-',/-0, . . .,* 



span the space, P/c, of polynomials of degree at most k and are linearly independent. Hence, if 
p(x) = do + . . .-{-OnX" (dj are arbitrary numbers) and k ^ n then 

p(x) = J &}*>(* )*J(1 -*)*"' (8) 

where 

« (!) 

"' is) 

with the assumption that a s = for s > n. 

We call (8) the generalized Bernstein form of p. 

THEOREM 2: //" p(x)= a + . . . + a n x n has real coefficients then for each k ^ n we have (putting 
K={0, . . ., k}) 

b(k)= nip ty k > ^ m ^ M ^ nip b^= B(k). 

PROOF: In view of (8) for each xel, p{x) is a convex combination of bb k \ . . ., b^K 
Following S. Bernstein, we can bound the discrepancies B(k)—M and m — b(k) as follows. 
We recall that if f(x) is defined on / 



*«>-M) (})*»-' 



j=o 
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Suppose ^ s ^ n. Then Bk(x s ; x) —x s ePk and hence 

B k (x s ;x)-x s =^8j(s) ( k )xJ(l-x) k -K 
j=o \J ' 

Theorem S: If k^n^l 



Sj( s ) ^ \ J J = °> ■ • •» k; s = 0, . . ., n. 



PROOF. Since Bk(l; x) = 1 and B/c (#; x) = x we have 

8 j (0)=8 i (l)=0,j=0, . . . 9 k. 

Henceforth, we assume s ^ 2. 
(i)0^y<s. 

In view of Lemma 1, we have 



^)=(f) s 



S-1V /5-l\ 2 (5-1 



(ii)2«s«jf 



8j(s)= f 



A« jUA-5)!_/jV K/-1) • • • U-(s-l)) 



k) (j-s)\k\ \k) k(k-l) . . . (k-(s-l)) 



\tji_tm 



1-1 



5-1 



Applying the mean value theorem to (1 — x) 8 ' 1 , we obtain 



1 _£Zi)-^i£Zil! i 



hence 



8j(s) 



(s-i) 2 /j\- % ^(s-iv 



-l) 2 (j\' 

k \kj 



Remark: If k > 2, (10) can be improved slightly to 

l\ (5-D 2 
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(10) 



(11) 



(12) 



since: 






















(i) If j 


<s, 
























8; (*) 


;< 


(s 


k 


1 
' k 


-(>■ 


-1)^ 


"I) 2 





Sfc(5)=0, 

and so (11) can be restricted to j < k leading to 

(5-1) 2 / k-iy (s-i) 2 ( k-\ 



3j{s) ^ k \ k ) k \ k 

THEOREM 4: //p(x) = a + . . . + a n x n , and k ^ n^ 1 then for j = 0, . . ., k 



satisfies 



p(^)-b< k) = 8 j (13) 



l«j|«£ (14) 



where 



//k^2, J/ierc 



A = £ (s-i) 2 Ia s |. (15) 

s=2 



|8,|<A^,j = 0, . . ., k. 



k ( k\ 
Proof: Since B k (p; x) —p{x) = ^ 8jl . j xi(l — x) k ~J, while in view of Theorem 3 

k0»; *) -pU) = 1 1 a*W ( *W(i -*)*-•>, 



#A: 



then 

k 
"o 



8 J = 2 «A(5), 



and the Theorem follows from (10) and (12). 
Corollary: IJ a , . . ., a n are real, then 



B(k)-M^|; m-b(k)^^. (16) 
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//k>2 



B(k)-M^A^;m-b(k)«A^- (17) 



Thus, B (k) converges to M as k-* °°, and to determine how large we need to choose k in order 
to be a given distance away from M, we need only consult (16) or (17) (analogously for b(k) and m). 
Cargo and Shisha (see footnote 1) give a difference table method for calculating b{ n \ . . . , b\p} x which 

can be used equally well to find Wf\ . . ., b { £\ (b^ ) = ao; b^ ) = ao-\- . . . + a n , always). 

Remark: Once a k has been chosen and the numbers 6 ( . fc) ordered according to size, it may 
be possible to improve on B(k) as an estimate for M, by means of a final correction. Suppose 

B(k) = b<*\ i^O, k, 

and 

max M*> = 6<*><6<*>. 



Put 



Pj(*) = (*)*>(l-x)*-i,./=0, 



LI1CI1 U ^= f /^ 


x ; \ 1 1UI clll A,tl , clIlU 

p(x)= Pi (x)b^+(l- 


- Pi {x)) ^bf( k 


/ 1-PiW 


Since 




] 


Pj(x) ^ Q 
-Pi(x) 




and 




2 


PjM _ 1 

1-Pi(x) 




we conclude 


from (18) that 










p{x)^ Pi (x)b^+(l- 


- Pi ( x ))b^=(b^- 


-ty*>)p,(*) +£<*>. 


Since Pi(x) 


assumes its maximum 


on / at 


x=ijk, we obtain. 


finally, 




M^ Pi ( 


i)»w+(i- w (i))^<fl(t). 


An analogous result holds for m. 
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(18) 



3. The Complex Case 

Suppose the coefficients a<>, . . ., a n are complex numbers. Let C be the convex hull of 
p(I) = {p(x) : *e/}, and let C k be the convex hull of {6 ( *>, . . ., &<*>}. Then in view of (8), we have, 
for each A: ^ ^ 

c c c*. 

We have 

LEMMA 2: For eac/i k ^ n 

Ck+i £ C k . 

PROOF: It is an easy consequence of (9) that b^ ) = ao, b^ ) = ao+ . . . +a„ and 

*5r i, -( 1 -*+i) l ? ,+ *+i%^- 1 *■ 



Thus 



bf+»eC k J=0, . . ., *+] 



and the lemma follows. 

Note that C and C*, fc = /i, ra+1, . . . are each compact sets in the plane, and that 



CQ (S C k . 

k=n 



(19) 



Indeed, 
Theorem 5: 



c= A c 



k- 



PROOF. Suppose there exists 



ze ff C/c 

k=n 



and zj^C. Since C is compact, dist (z, C) = d > 0. Choose k ^ n so large that 04//c) < d. Now, 
there exist \o, • • ., ^k satisfying kj ^ 0, SXj = 1 such that 



z=£ Mf , 



j=o 



hence, in view of (13) and (14) 



and 



k / j\ k 

j=o Kh/ j=o 



tMi) 



<d. 



(20) 
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But 



± *({)«. 



j=0 

and (20) contradicts dist (z, C) =d, thus proving the theorem. 
Note that the same proof shows that 

dist (C*,C)<~ 
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